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ABSTRACT 



The spectrometer SPI on INTEGRAL allows for the 
first time simultaneous imaging of diffuse and point-like 
emission in the hard X-ray and soft gamma-ray regime. 
To fully exploit the capabilities of the instrument, we im- 
plemented the MREM image deconvolution algorithm, 
initially developed for COMPTEL data analysis, to SPI 
data analysis. We present the performances of the algo- 
rithm by means of simulations and apply it to data accu- 
mulated during the first 2 mission years of INTEGRAL. 
Skymaps are presented for the 1809 keV gamma-ray line, 
attributed to the radioactive decay of 26 A1, and for contin- 
uum energy bands, covering the range 20 keV - 3 MeV. 
The 1809 keV map indicates that emission is clearly de- 
tected by SPI from the inner Galactic radian and from the 
Cygnus region. The continuum maps reveal the transi- 
tion between a point-source dominated hard X-ray sky to 
a diffuse emission dominated soft gamma-ray sky. From 
the skymaps, we extract a Galactic ridge emission spec- 
trum that matches well SPI results obtained by model fit- 
ting 0,0. By comparing our spectrum with the cumu- 
lative flux measured by IBIS from point sources, we find 
indications for the existence of an unresolved or diffuse 
emission component above ^100 keV. 

Key words: Image deconvolution; 1809 keV line emis- 
sion; Galactic continnum emission. 



1. INTRODUCTION 



The image deconvolution of gamma-ray data is con- 
fronted with several problems that are specific to the 
domain. First, the data are dominated by a strong in- 
strumental background, generated by the interaction of 
cosmic-rays with the telescope detectors and the sur- 
rounding satellite materials. For SPI, typical signal-to- 
background ratios are of the order of ~1%. Hence, the 
deconvolution seeks to explain only a tiny fraction of the 
registered events, while the bulk of the data have to be ex- 



plained by a model of the instrumental background char- 
acteristics. Second, the information measured by the tele- 
scope for each registered photon is only indirectly related 
to its arrival direction. For SPI, the general data space is 
3 -dimensional, spanned by the pointing number, the de- 
tector identifier, and the measured energy. Analysing the 
data in a single energy band reduces the date space to two 
dimensions, which differ however considerably in nature 
from the image space which is spanned by longitude and 
latitude. Third, measured count rates are generally very 
low, leading to non Gaussian distributions for the statisti- 
cal measurement uncertainties. The consideration of the 
Poissonian nature of the data is in many cases mandatory 
to derive meaningful error estimates. 

All these particularities explain why imaging the gamma- 
ray sky with SPI aboard INTEGRAL is difficult, and 
why classical image reconstruction methods rapidly show 
their limitations in this domain. Statistical noise easily 
propagates into the reconstructed sky intensity distribu- 
tions, leading to artefacts that can severely mislead the 
observer. One may seek to suppress this statistical noise 
using regularisation methods, premature stopping of the 
iterations, and/or image smoothing 
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tricks go however in general at the expense of flux sup- 
pression and a loss in the angular resolution. 

We discuss in this paper the application of an alterna- 
tive method for image deconvolution of INTEGRAL/SPI 
data which is based on a multiresolution analysis of the 
data using wavelets. The method is an implementation of 
the Multiresolution Expectation Maximisation (MREM) 
algorithm that has been initially developed to perform 
image deconvolution of gamma-ray data that where ac- 
quired by the COMPTEL telescope aboard the Comp- 
ton Gamma-Ray Observatory (CGRO) [9]. The algo- 
rithm is implemented as the program spi_obs_mrem 
(version 3.4.0) and is publicly available at the web site 
|http : / / www .cesr.fr/ ~ jurgen/ isdc| We will 
demonstrate that MREM is equally well suited for imag- 
ing of point source and diffuse emission, that the algo- 
rithm is able to cope with large exposure gradients over 
the imaged field, and that it allows at the same time reli- 
able estimates of gamma-ray intensities from the result- 



ing images. MREM is a convergent algorithm, and it will 
be shown that its only sensitive parameter is the wavelet 
thresholding level which has a well defined scaling, al- 
lowing for an objective a priori choice of its value. 



2. DECONVOLUTION ALGORITHMS 



2.1. Richardson-Lucy algorithm 

The starting point of MREM is the Richardson-Lucy 
(RL) algorithm dH 

which is widely used for the de- 
convolution of astronomical images, and which, in partic- 
ular, has been also successfully employed for the analysis 
of gamma-ray data of CGRO ISJ O- Starting from an 
initial estimate /? for the image which we usually take 
as a grey image of negligible intensity^], RL iteratively 
improves this estimate using the relation 

f k+i = f k + 6 k (1) 

with 
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being the image increment of iteration k and 
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being the predicted number of counts in data space bin i 
after iteration k, where 

• i is the index of the data space, running from 1 to N, 

• j is the index of the image space, running from 1 to 

M, 

• k counts the number of iterations, starting with for 
the initial image estimate, 

• f!f if the intensity in image pixel j after iteration k 
(units: ph cm _2 s~ 1 sr _1 ), 

• Rij is the instrumental response matrix linking the 
data space to the image space (units: counts ph -1 
cm 2 s sr), 

• m is the number of counts measured in data space 
bin i (units: counts), 

• hi is the predicted number of instrumental back- 
ground counts for data space bin i (units: counts). 

Convergence of this algorithm may eventually be very 
slow, in particular if the instrumental background, &j, 
is the dominating term in Eq. (01. Further, the model 

1 This choice corresponds to having no a priori assumption about 
the possible spatial distribution of the emission. 



of the instrumental background may be presented in a 
parametrised form, e.g. ..,(3 m ), where the back- 

ground model parameters f3% , .., (3 m are to be determined 
along with the image deconvolution process. 

Convergence acceleration and simultaneous background 
model parameter fitting is achieved by replacing Eq. (TJ 
by 

f- +1 = // + \ k &l (4) 

and by optimising the acceleration factor A fe together with 
the background model parameters fix,.., j3 m by maximis- 
ing the log-likelihood function 

N 
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with 
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To conserve positivity of the image pixels during the ac- 
celeration process (a property which is inherent to the 
standard RL algorithm) we require X k Sj > —ft for all 
image pixels j. 



2.2. MREM algorithm 

The RL algorithm as implemeted by Eqs. |U |2] E] and 
[6] does not explicitely take into account the correlation 
of intensities in neighbouring image pixels. In the case 
of weak signals, individual pixels are therefore poorly 
constrained by the data and statistical fluctuations easily 
propagate in the reconstruction and quickly amplify. One 
approach to limit this noise amplification is the applica- 
tion of a smoothing kernel to the image increments Sj, 
i.e. to replace 

M 
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As an example, the 511 keV allsky map presented in ifTTll 
has been obtained by the RL algorithm using a 5°x5° 
boxcar average. Together with premature stopping of the 
iterations, this resulted in a basically noise free image of 
the Galactic 511 ke V intensity distribution. 

The drawback of this approach is (1) the effective lim- 
ited angular resolution imposed by the smoothing proce- 
dure and (2) the arbitrary stopping of the iterations. In 
the MREM approach, we circumvent these problems by 
replacing Eq. |7]by the wavelet filtering step 

M' 
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where 
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are the wavelet coefficients of the image increment, W sr 
is a M' x M matrix that present the wavelet transform 
(and W' 1 the inverse transform), and 

,; = {;: i<* 

is a threshold vector for the wavelet coefficients c k . The 
threshold vector represents the essence of the filtering 
step in that it zeros all those wavelet coefficients that 
fall below a given threshold t k . The idea behind this 
kind of filtering is that significant image structures are 
represented by large wavelet coefficients while insignif- 
icant structures are represented by small coefficients. 
The wavelet transform operates globally on the image 
and therefore introduces pixel-to-pixel correlations on all 
possible scales, and the careful selection of t k allows us 
to keep only structures for which significant evidence is 
present in the data. Once all significant structures have 
been extracted, the algorithm stops. 

We estimate t k s for each iteration using Monte-Carlo sim- 
ulations. For this purpose we create simulated image in- 
crements by replacing in Eq. [2] the measured number of 
counts, Hi, by the mock data n; L which have been drawn 
from the current estimates e k using a Poisson random 
number generator. Since rij are a possible realisation of 
e k , the resulting image increment reflects now the impact 
of statistical noise fluctuations only. Transforming this 
fake image increment in the wavelet domain using Eq. [9] 
gives then possible values for all wavelet coefficients in 
the presence of pure statistical noise. Repeating this pro- 
cess a sufficient number of times allows then to estimate 
the root-mean-square, a k , for each coefficient. We use 
then these estimates and set the thresholds to 

t k s =axa k s (11) 

where a may be interpreted as the number of standard 
deviations of noise that are allowed during the recon- 
struction. We found satisfactory results for a «3.5. 100 
Monte Carlo simulations per iteration were found suffi- 
cient for reliably estimating the noise levels a k . 



3. SIMULATIONS 



To illustrate the performances of MREM we present two 
simulated reconstructions in Figs. [Hand [2] Both simu- 
lations were based on the SPI pointing sequence of the 
first 2 years of observations, spanning the orbital revolu- 
tions 19-269. In both cases, we adopt a realistic model 
for the instrumental background that was based on the 
measured data. To simulate the gamma-ray sky, we con- 
volve model intensity distributions with the instrumental 
response, and add the resulting data space model to that 
of the instrumental background. Mock data are then gen- 
erated from these models using our tool spi_obs_sim, 
which employs a Poisson random number generator to 
create a realisation of measured events that is statistically 
compatible with the model. 



Our first example illustrates the MREM reconstruction of 
point source emission (Fig.[T]i. As sky model we used the 
IBIS catalogue of 49 sources found in the 100-150 keV 
energy band [2]. The mock dataset has been created and 
the MREM reconstruction has been performed for the 
same energy band, but now for the SPI instrument which 
is less sensitive with respect to IBIS. Consequently, from 
the 49 simulated sources MREM reaveals only the 15 
most brightest ones in the reconstruction. In addition, the 
modest angular resolution of SPI of ^3° leads to some 
source confusion towards the densely populated Galac- 
tic centre region. Besides the 15 point sources, no other 
emission features are present, and in particular, the map 
is free from image noise. 

Our second example illustrates the MREM reconstruction 
of diffuse emission (Fig. [2]). As sky model we've chosen 
the Dirbe 240 /im cold dust emission map which has been 
shown to provide a reasonable representation of the spa- 
tial distribution of 1809 keV gamma-ray line emission at- 
tributed to the radioactive decay of 26 A1 [ 10]. In order to 
highlight details of the reconstruction process, we scaled 
the intensity of the skymap to approximately 5 times of 
that expected from 1809 keV line emission. The simula- 
tion and reconstruction was performed for a 7 keV wide 
energy band centred on 1808.5 keV. MREM clearly picks 
up the basic features of the input map: bright and narrow 
Galactic ridge emission plus several localised emission 
features in the Cygnus, Carina and Vela regions. Low- 
intensity emission that is coded in reddish color in the 
model map falls apparently below the SPI detection limit, 
and consequently is clipped by the wavelet filter during 
the deconvolution. Also for this case, the reconstruction 
appears to be free of image noise. 

Another interesting feature of the MREM reconstruc- 
tion is the sharpening of the Galactic ridge emission to- 
wards the regions of high-intensity in the inner Galaxy. 
Since the signal in this area is detected at higher statis- 
tical significance, the data allow to constrain the inten- 
sity profile more precisely, and the reconstructed profile 
becomes sharper. Further away from the Galactic cen- 
tre, say beyond ±40° Galactic longitude, the weaker sig- 
nal intensity makes the data less constraining, and conse- 
quently, the reconstructed profile becomes broader. This 
behaviour is well known in statistics as the bias-variance 
trade-off, which generally occurs when we try to separate 
signal from noise using smoothing or correlation tech- 
niques. 



4. 1809 KEV GAMMA-RAY LINE EMISSION 



We now apply the MREM algorithm to real data. The 
dataset we used is the same as that employed for the simu- 
lations: the first 2 years of INTEGRAL/SPI observations, 
spanning the orbital revolutions 19-269. As usually in 
our analyses, we screened the data for solar flares, excep- 
tionally high countrates, and other possible problems. 

Our first goal is to image the sky in the 1809 keV gamma- 
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Figure 1. Model ( top) and MREM reconstruction (bottom) of a mock dataset of point sources. The simulated sky model 
is a compilation of IBIS point-sources detected in the 100-150 keV band /HI/. For the purpose of better visibility, the 
emission for each source in the model has been spread over a diameter of 2 degrees around the source position. Image 
reconstruction has been performed using the SPI response in the 100-150 keV energy band. 



ray line attributed to the radioactive decay of 26 A1. Simi- 
lar to the simulation, we analyse the data in a 7 keV wide 
energy band centred on 1808.5 keV. To fully exploit the 
sensitivity of SPI we use both single and double events. 
The resulting MREM image is shown in the top panel of 
Fig. [5] For comparison, we also show the results of the 
RL reconstruction after iterations 30 and 60 in the mid- 
dle and bottom panels of Fig. [3] To reduce image noise 
in the RL reconstruction we applied a boxcar average of 
15°xl5° to smooth the image increment after each iter- 
ation (cf. Eq. 0. The images may also be compared to 
the RL image presented in [6] which uses an exponential 
disk as initial image estimate /? and for which no boxcar 
averaging has been performed. 

The MREM image indicates that 1809 keV line emis- 
sion is clearly detected by SPI from the inner Galactic 



radian (say between ±45° Galactic longitude) and from 
the Cygnus region (situated in the Galactic plane at ~80° 
longitude). The inner Galaxy emission shows some lo- 
cal maxima, but those are probably of low statistical sig- 
nificance. Indeed, when we increase the wavelet thresh- 
old level from our standard value of a=3.5 to 4 most 
of these maxima disappear, suggesting that the features 
were probably linked to statistical noise fluctuations. 

COMPTEL observations have indicated 1809 keV line 
emission also from the Carina and Vela regions, located 
along the Galactic plane at longitudes ^286° and ^266°, 
respectively |8|,|5|]. The RL reconstruction shows indeed 
a hint for emission from Carina (visible after iteration 
60), but the large amount of structure in the image to- 
wards the inner Galaxy - which is due to statistical noise 
fluctuations - indicates that the significance of the Carina 
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Figure 2. Model ( top) and MREM reconstruction (bottom) of a mock dataset of diffuse emission. The simulated sky model 
is the Dirbe 240 pm cold dust emission map, scaled to an intensity that corresponds to about 5 times the intensity of the 
Galactic 1809 keV line emission. Image reconstruction has been performed using the SPI response in a 7 keV wide energy 
band centred on 1808.5 keV. 



feature is weak in our analysed SPI dataset. So far our 
imaging analysis shows no indication for 1809 keV line 
emission from the Vela region, but also in the COMPTEL 
data, the statistical significance of this feature was among 
the weakest (see also Schanne et al., these proceedings). 



age similiar to that shown after iteration 60. In fact, at this 
point the image increment smoothing did prevent any fur- 
ther improvement of the likelihood, so the iterations were 
aborted. 



It is interesting to note that the RL image after iteration 30 
appears similar to the MREM image. The heavy boxcar 
smoothing ofl5°xl5° that we applied during the RL de- 
convolution plays apparently a similar role as the wavelet 
filter in that it suppresses high-frequency image noise 
during the reconstruction. However, the RL iterations do 
not stop at this point, and even such heavy smoothing 
does not prevent the growth of image noise in the recon- 
structions with proceeding iterations. Although we have 
not formulated an explicit stopping criterion for RL, the 
algorithm did indeed stop after 85 iterations with an im- 



5. GALACTIC CONTINUUM EMISSION 



Our second goal is the imaging of the hard X-ray and 
soft gamma-ray emission in broad continuum energy 
bands. This energy range is a transition region where 
at least three emission components are present: (1) soft 
(^300 keV) point source emission from Galactic X-ray 
binaries and AGNs, (2) diffuse bulge dominated emis- 
sion up to 511 keV attributed to ortho-positronium an- 
nihilation, and (3) hard (Jl300 keV) extended Galactic 




Figure 3. Comparison of MREM (top panel) and RL (middle and bottom panel) reconstructions of 1809 keV data. For 
RL we show the result after iteration 30 and 60. To reduce the image noise, we applied a boxcar average ofl5°xl5°for 
the RL reconstruction. 




Figure 4. MREM images of the hard X-ray sky in the continuum energy bands 20-25, 25-50, and 50-100 keV. 




Figure 5. MREM images of the hard X-ray and soft gamma-ray sky in the continuum energy bands 100-200, 200-300, 
and 300-500 keV. 
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Figure 6. MREM images of the gamma-ray sky in the continuum energy bands 514-1000 keV, and 1-3 MeV. 



plane emission of yet unknown nature, possibly linked 
to the interaction of cosmic -rays with the interstellar 
medium. The spectral characteristics of these 3 com- 
ponents have been already extensively studied with SPI 
HEU], but images of this transition region, in particular 
above <;100 keV have so far only been published for the 
ortho-positronium annihilation component B20I1 . 

In Figs.[4]-[6]we present the MREM reconstructions for 8 
energy bands situated between 20 keV and 3 MeV. The 
hard X-ray images below 100 keV nicely show the popu- 
lation of Galactic point sources, mainly X-ray binaries, 
of which we detect ^60 in the images below 50 keV 
and ^25 in the 50-100 keV image. No diffuse emis- 
sion is visible in the reconstructions, indicating that, if 
present, such emission should be at a rather low inten- 
sity level. In the 100-200 keV image we still count ^18 
point sources in the image, without a clear hint for diffuse 
emission. In the 200-300 keV image, a diffuse and struc- 
tured emission band along the Galactic plane replaces the 



region where the point sources were located at lower en- 
ergies. Emission maxima coincide well with the posi- 
tions of point sources found in the 100-200 keV image, 
and we therefore suggest that at least part of the 200- 
300 keV emission originates from faint and unresolved 
point sources. We cannot exclude, however, that a sub- 
stantial fraction of the emission in this energy band is in- 
deed of diffuse nature. 

The morphology of the emission changes drastically 
above 300 keV. The 300-500 keV image is dominated 
by Galactic bulge emission, which is very similar in mor- 
phology to the emission that is observed in the 511 keV 
gamma-ray line (see Weidenspointner et al., these pro- 
ceedings). The only point sources that are still clearly 
present at these energies are the Crab (near the anticen- 
tre) and the X-ray binary Cyg X-l. The bulge emission 
seems to show some flattening which may be taken as 
an indication of another, yet weaker, underlying emission 
component that follows the Galactic plane. 
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Figure 7. Comparison of the MREM inner Galactic radian spectrum (filled dots) to SPI model fitting results (crosses, 
[19]) and to spectra obtained by other instruments. The filled squares show the summed point source contributions from 
the second IBIS catalogue normalised to the same Crab spectrum as used for the SPI data. The open dots show the 

residual MREM flux after subtraction of the IBIS point source contribution. The open squares show the diffuse hard X-ray 
emission found by [12] using IBIS/ISGRI as a light-bucket. The triangles show the result obtained from COMPTEL data 
U7\l , the diamonds show the OSSE results ^/. The black line is the most recent GALPROP model prediction for galactic 
diffuse emission (Strong, private communication). 



Above 514 keV (i.e. at energies situated above the 511 
keV positron-electron annihilation line), the gamma-ray 
sky becomes rather dark (cf. Fig.|6]l. Only 2 point sources 
are visible in the 514-1000 keV and 1-3 MeV images, 
the Crab and Cyg X-l, while no emission can be seen 
from the Galactic plane. If we would plot the images to 
very low intensities we would start to see some irregular 
emission correlated with the Galactic disk, but MREM 
did not manage to extract this component clearly from the 
data. In other words, the extended Galactic disk emission 
that we know from spectral analysis to be present above 
514 keV in our SPI observations is too weak to be 

reliably imaged with the present dataset. 

However, we can determine the global intensity of the 
Galaxy at these energies by integrating the flux that is 
present in the MREM images over large regions. We 
show the result of such an integration in Fig. [7] In or- 
der to determine the global flux from the inner Galactic 
radian, we have chosen an integration region of ±30° in 
Galactic longitudes and ±10° in Galactic latitudes. We 
normalised the spectrum to that of the Crab by assuming 
for the latter a broken power law with a break energy of 
1 17 keV, power law indices of -2.04 and -2.37 for below 
and above the break, respectively, and a normalisation of 
10.78 at 1 keV. We derived this prescription by fitting SPI 
data of the Crab over the 25-500 keV energy band. 

For comparison we show in Fig.|7]also the spectrum that 



we obtained from SPI data using a model fitting proce- 
dure IllSfll . normalised to the above mentioned Crab spec- 
trum. The MREM result matches quite nicely the spec- 
trum obtained by model fitting. The advantage over the 
model fitting approach is that the MREM results are un- 
biased towards assumptions about the spatial distribution 
of the emission, and this difference may eventually ex- 
plain the slightly higher fluxes obtained by MREM in the 
— 100 - 200 keV range. The disadvantage is that MREM 
does not provide an estimate of the statistical uncertain- 
ties related to the flux measurements. Hence, both meth- 
ods are best used in conjunction, providing thus reliable 
estimates of both the spatial and spectral distribution of 
the emission. 

We also show in Fig. [7] the spectrum that is obtained by 
summing the flux of all point sources listed in the sec- 
ond IBIS catalogues [3, 2] within the integration region. 
We again use the Crab normalisation as for the MREM 
points to convert the catalogue fluxes (quoted in mCrab) 
into physical units. For energies below 100 keV, the 
MREM fluxes match fairly well the point source fluxes 
determined with IBIS/ISGRI, indicating that most of the 
emisison in this band can indeed be attributed to point 
sources 111311 . Subtracting the IBIS/ISGRI point source 
flux from the SPI MREM measurements provides an esti- 
mate of the unresolved (or diffuse) emission component, 
which is shown as open dots in Fig. [7] Below ~ 40 keV, 



this estimate is consistent with the findings of |[T1 who 
report the detection of an unresolved emission compo- 
nent using IBIS/ISGRI that they attribute to a large pop- 
ulation of weak CVs (open squares). Above ~ 40 keV, 
our measurements indicate an additional component that 
is not resolved by IBIS/ISGRI into point sources. While 
above ~ 300 keV ortho-positronium annihilation towards 
the galactic bulge may explain partly this discrepancy, the 
excess emission between ~ 40 — 300 keV can only be at- 
tributed to an additional unresolved or diffuse emission 
component. 

We also show in Fig. [7] spectral points obtained by the 
COMPTEL [17] and OSSE H telescopes, but we did 
not attempt a cross-calibration using the Crab for these 
instruments, which may partially explain some apparent 
discrepancies (even between both instruments). Above 
511 keV, the galactic emission is believed to originate 
mainly from cosmic-rays interacting with the interstel- 
lar medium and the interstellar radiation field, and we 
show in Fig. [7] the latest prediction of the GALPROP 
cosmic-ray propagation model as solid line (Strong, pri- 
vate communication). Apparently, GALPROP explains 
relatively well the spectral shape of the unresolved (or 
diffuse) emission component, both below and above the 
511 keV annihilation line, yet the flux level is underesti- 
mated by a factor ^2 — 3. Whether this underestimation 
is due to a yet unrecognised diffuse emission mechanism 
or due to a yet unresolved weak source population re- 
mains to be seen. 

Fig. [7] illustrates that MREM images may be exploited 
quantitatively to derive source spectra. However, such an 
analysis can not replace a detailed model fitting approach 
which allows to assess the significance of emission com- 
ponents and to determine the statistical uncertainties in 
the flux estimates. MREM spectra should therefore be 
understood as complementary to the model fitting ap- 
proach in that they do not assume a specific spatial dis- 
tribution for the emission components. MREM images 
should hint towards the components that are needed to 
satisfactorily describe the data, model fitting should then 
be used to assess how significant these components in- 
deed are. 



6. CONCLUSIONS 



The MREM image reconstruction algorithm presents an 
interesting alternative for the deconvolution of SPI hard 
X-ray and gamma-ray data. The main features of this 
method are (1) the suppression of image noise in the re- 
constructions and (2) the convergence towards a well- 
defined solution. The application of MREM to simu- 
lations nicely illustrates that these properties are indeed 
reached by the present implementation. 

We employed the algorithm to present SPI images in 
the 1809 keV gamma-ray line from 26 A1 and in contin- 
uum bands covering energies from 20 keV up to 3 MeV. 
The 1809 keV image confirms the basic emission fea- 



tures that have been established by COMPTEL, namely 
the strong Galactic ridge emission from the inner Galaxy, 
and a prominent emission feature in Cygnus. The con- 
tinuum images reveal nicely the transition from a point 
source dominated hard X-ray sky to a diffuse emission 
dominated soft gamma-ray sky. Imaging the sky above 
511 keV is still a challenge, even for MREM, due to 
the low intensity of the Galactic emission and the corre- 
spondingly small signal that is present in the data. Nev- 
ertheless, MREM allowed to derive spectral points for 
the Galactic emission up to 1 MeV, which are in fairly 
good agreement with previous measurements by OSSE 
and COMPTEL. 

The accumulation of data beyond the 2 years database 
used for this analysis will enable more detailed imaging 
of hard X-ray and soft gamma-ray emission throughout 
the spectrum that is accessible to SPI. We expect to im- 
age the diffuse Galactic emission that is underlying the 
strong ortho-positronium continuum in the intermediate 
300-511 keV band, and we hope to visualise also the 
spatial distribution of Galactic continuum emission above 
511 keV. Imaging those emissions will be crucial for the 
determination of their origin, which theoretically, is still 
poorly understood tl% . 
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